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Iron-chromium alloys, the base components of various stainless steel grades, have numerous tech¬ 
nologically and scientifically interesting properties. However, these features are not yet sufficiently 
understood to allow their full exploitation in technological applications. In this work, we investigate 
segregation, precipitation, and phase separation in Fe-Cr systems analysing the physical mechanisms 
behind the observed phenomena. To get a comprehensive picture of Fe-Cr alloys as a function of 
composition, temperature and time the present investigation combines Monte Carlo simulations us¬ 
ing semiempirical interatomic potential, first-principles total energy calculations, and experimental 
spectroscopy. In order to obtain a general picture of the relation of the atomic interactions and 
properties of Fe-Cr alloys in bulk, surface, and interface regions several complementary methods has 
to be used. Using Exact Muffin-Tin Orbitals method the effective chemical potential as a function 
of Cr content (0-15 at.% Cr) is calculated for a surface, second atomic layer and bulk. At ~10 at.% 

Cr in the alloy the reversal of the driving force of a Cr atom to occupy either bulk or surface sites 
is obtained. The Cr containing surfaces are expected when the Cr content exceeds ~10 at.%. The 
second atomic layer forms about 0.3 eV barrier for the migration of Cr atoms between bulk and 
surface atomic layer. To get information on Fe-Cr in larger scales we use semiempirical methods. 

Using combined Monte Carlo molecular dynamics simulations, based on semiempirical potential, 
the precipitation of Cr into isolated pockets in bulk Fe-Cr and the upper limit of the solubility of 
Cr into Fe layers in Fe/Cr layer system is studied. The theoretical predictions are tested using spec¬ 
troscopic measurements. Hard X-ray photoelectron spectroscopy and Auger electron spectroscopy 
investigations were carried out to explore Cr segregation and precipitation in Fe/Cr double layer 
and Fe0.95Cr0.05 and Feo. 85 Cro .15 alloys. Initial oxidation of Fe-Cr was investigated experimentally 
at 10~® Torr pressure of the spectrometers showing intense Cr 2 03 signal. Cr segregation and the 
formation of Cr rich precipitates were traced by analysing the experimental spectral intensities with 
respect to annealing time, Cr content, and kinetic energy of the exited electron. 

PACS numbers: 68.35.bd, 68.35.Dv, 68.47.De, 71.15.Nc 


I. INTRODUCTION 

Iron-chromium alloys have many technologically im¬ 
portant and scientifically interesting properties.^ On the 
other hand, growing technological challenges are faced in 
designing multifunctional steels. Developing of advanced 


steels to meet different standards simultaneously, e.g. 
high strength, proper workability and ductility, excellent 
corrosion resistivity, and specific magnetic properties all 
in various ambient conditions, require breakthrough in¬ 
novations and cutting edge research.^ Iron and chromium 
are typical examples of a ferromagnet and an antiferro- 
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magnet, respectively. As a function of composition and 
structure the magnetic properties of Fe-Cr vary consid¬ 
erably, e.g. spin glass^ and giant magnetoresistance^ fea¬ 
tures are found in these systems. The crystallographic 
properties of Fe-Cr are also peculiar. Although Fe-Cr 
has a body centered cubic (bcc) based structure within 
the whole composition range, there exist both stable 
and metastable composition regions in the phase dia¬ 
gram. The experimental phase diagram of Fe-Cr at 300 
°C shows the miscibility gap beginning from 5-10 at.% 
and extending to 90-95 at.% Cr. Within this miscibility 
gap two domains exist, regions of spinodal decomposition 
and nucleation and growth. The latter is located at the 
outskirts of the miscibility gap extending about 20 at.% 
on each side.® On aging Fe-Cr often undergoes trans¬ 
formation to either high-temperature cr-phase or separa¬ 
tion into Fe-rich (a) and Cr-rich (a’) phases.In 20 
at.% Cr alloy at 773 K temperature after 50 h anneal¬ 
ing precipitates have been observed to occupy 2 % of the 
alloy volume.® The microscale changes in the crystallo¬ 
graphic properties may induce considerable changes at 
macroscopic level, e.g. the ’475 °C embrittlement’^^’^^, 
has significant effect on the mechanical properties of cer¬ 
tain steel grades. Fe-Cr is the base component in many 
stainless steel grades due to the beneficial properties of 
chromium. Certain amount of chromium makes an iron 
alloy corrosion resistant^^. At ambient conditions a thin 
and transparent film of chromium oxide rapidly forms 
on the open surface of the alloy preventing further oxi¬ 
dation and blocking corrosion. The corrosion resistance 
of the ferritic stainless steels increases abruptly by sev¬ 
eral orders of magnitude when the Cr content in bulk 
reaches ^10 at.% level.This oxidation-related experi¬ 
mental threshold of Cr content in bulk coincides with the 
calculated reversal point of the relative magnitudes of the 
Fe and Cr chemical potentials in bulk and surface of the 
Fe-Cr alloys. This reversal of the relative chemical po¬ 
tentials enables the outburst of Cr on the otherwise pure 
Fe surface found exclusively in the case of low-Cr Fe-Cr 
alloys.Therefore, at ambient conditions Cr 2 03 is easily 
formed on the surface of Fei_xCrx (x > 0.1) alloys. Fur¬ 
thermore, due to the strong tendency of Cr to segregate 
to the Fe-Cr/Cr 203 interface,^® there is an additional 
driving force for a Cr 203 layer to grow until the surface 
oxide reaches the protective nanometer scale thickness 
preventing the further oxidation of the material. Due 
to their technological importance and challenging open 
questions the surfaces of Fe-Cr alloys have recently re¬ 
ceived considerable scientific attention.^^44,16-21 

In this work, we investigate the physical conditions and 
possible realizations of segregation, precipitation, and 
phase separation in Fe-Cr systems. To get a comprehen¬ 
sive picture of the state of Fe-Cr as a function of com¬ 
position, temperature and time the present investigation 
combines Monte Carlo simulations using semiempirical 
interatomic potential, ab initio total energy calculations, 
and experimental spectroscopy. Using several comple¬ 
mentary methods it is possible to get a more reliable 


picture of the interactions between Fe and Cr atoms and 
explain the consequences of the atomic interactions for 
the properties of Fe-Cr alloys. Ah initio methods are 
used to get the atomic scale energetics as a function of 
the concentration of the alloys, thermodynamics of large 
scale systems is obtained by Monte Carlo method, and 
experimental spectroscopy is used to probe the concen¬ 
trations and atomic structure of real Fe-Cr systems as a 
function of Cr content and annealing time. However, one 
should carefully take into account the scope of validity of 
each method. Keeping this in mind our aim is to give an 
extensive picture of the atomic structure in bulk, surface, 
and interface regions in Fe-Cr systems. The rest of the 
paper is divided into two main sections and conclusions. 
The research methods are briefly reviewed in Section II 
and the results are presented and discussed in Section 
III. 


II. METHODS 

A. First-priciples calculations 

The ab initio total energy calculations are based on 
the density functional theory^^^a were performed us¬ 
ing the Exact MufHn-Tin Orbitals (EMTO) method.^'^^5 
The basis set includes s, p, d, and / orbitals. The gen¬ 
eralized gradient approximation in the PBE form was 
used for the exchange-correlation functional.^® The to¬ 
tal energy was calculated using the full charge-density 
technique.^®4’^ The alloys were simulated as substitution- 
ally disordered ferromagnetic bcc phase using the coher¬ 
ent potential approximation (CPA) which provides the 
continuous scanning of the concentration of the alloy.^® 
The calculated equilibrium lattice constant was used for 
each composition. The EMTO approach in combination 
with the CPA has been applied successfully in the theo¬ 
retical study of various structural and electronic proper¬ 
ties of alloys and compounds^® demonstrating the level 
of accuracy and efficiency needed also in the present in¬ 
vestigation. Eor more details of the electronic structure 
calculations we refer to our earlier work. ^'^49 

The basic quantities used in the present study are sur¬ 
face energy (Agurf), segregation energy (Agegr), chemical 
potential (p) and mixing enthalpy per atom (if„iix)- The 
effective chemical potential (A/i®) and the slope of the 
mixing enthalpy of bulk Fe-Cr are related within a simple 
relation,®® 

A/i® = (/rpe - «-+ constant, (1) 

where x is the atomic fraction of Cr (Acr/(AFe + Acr), 
Ape and Acr are the number of Fe and Cr atoms in the 
investigated system, respectively). The surface energy is 
dehned as the energy needed to form a new surface per 
the formed new surface area. 

Agiab Abulk /r,\ 
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where -©slab is the energy of the slab system with two 
surfaces, both having the area A. ifbuik is the energy 
of the bulk system having the same amount of atoms as 
the slab system. The segregation energy of Cr from a 
region A to a region B is defined as the energy needed to 
transfer a Cr atom from A to B and an Fe atom from B 
to A. 

(3) 


B. Large-scale Monte Carlo simulations 

Because ab initio simulations for large systems and for 
longer time scales are not possible we performed Monte 
Carlo simulations for Fei-xCrj, alloys to test the ab initio 
predictions and to investigate the thermodynamic ground 
state of the alloys in larger scales. 

Due to the time consuming simulations and large sys¬ 
tems, ab initio methods can not be used. Thus the in¬ 
teratomic interaction was modeled by a semiempirical 
potential, namely the two-band embedded atom model 
(2BEAM) which is designed to reproduce the mixing en¬ 
thalpy of the Fe-Cr alloy.^^ One should note that the 
2BEAM potential is optimized mainly for bulk proper¬ 
ties. 

Precipitation and segregation of chromium in finite 
temperatures and in large systems was studied using 
the Monte Carlo (MC) method where possible moves in¬ 
cluded atom displacements and exchange of types (Ee or 
Cr) of a pair of atoms. Displacements were performed 
with short sequences of molecular dynamics (MD) sim¬ 
ulations in the canonical ensemble. Using MD was ob¬ 
served to be more efficient in moving atoms than the 
conventional Metropolis algorithm with atomic displace¬ 
ments. All Monte Carlo-molecular dynamics (MCMD) 
simulations were performed in NVT ensemble with the 
proper value of the lattice constant obtained from sepa¬ 
rate NPT simulations as a function of temperature and 
chromium concentration. It should be emphasized that 
the MCMD calculations are pure equilibrium simula¬ 
tions; there is no kinetics involved. 

The MCMD method was used to study the near sur¬ 
face structure of the Fe-Cr alloy and the structure of a 
iron-chromium interface in a layered system. In the sim¬ 
ulations a system of size 86 x 86 x 86 with 54000 atoms 
was used. The lengths of the simulations varied from 
80000 to 120000 MC steps and results were calculated 
by taking the averages of roughly 40000 last simulation 
steps. For the surface studies boundary conditions were 
applied in the x and y directions while leaving the two 
(001) 2 surfaces open. In the case of layer structure sim¬ 
ulations periodic boundary conditions were applied in all 
three directions. The interface orientation was (001). 

Surface structure simulations were performed in tem¬ 
peratures of 300, 500, and 700 K and interface simula¬ 
tions in 300 and 700 K. 


C. Experiments 

The Fe/Cr bilayer was grown by electron beam physi¬ 
cal vapour deposition from elemental Fe and Cr on a Si 
substrate. The top Fe layer of the Fe/Cr/Si sandwich was 
grown to about 50 nm thickness to protect the sample 
from contamination and mechanical failure. The thick¬ 
ness of the Fe film was checked by sputtering. The Fe-Cr 
alloy samples were prepared by induction melting under 
argon flow from elemental components. The Cr bulk con¬ 
centrations are 5 and 15 at. % for samples Fe0.95Cr0.05 
and Feo. 85 Cro. 15 , respectively. The concentrations 5 and 
15 at. % Cr were selected to encompass the interesting 
concentration region of the onset of the corrosion resis¬ 
tance (9-13 at. % Cr) in ferritic stainless steels and the 
ab initio prediction of the onset of the surface segregation 
of Cr in Fe-Cr alloys (8-9 at. % Cr).^"^ Elemental compo¬ 
nents of purity better than 99.99 % were used for all the 
samples. To start the investigations with fresh and un¬ 
oxidized samples sputtering and annealing were used for 
all the samples. The Fe layers of Fe/Cr/Si were cleaned 
in the UHV of analyzer chamber of spectrometer by an¬ 
nealing at 150 °C followed by 20 min argon sputtering. 
To clean the alloy samples Ar sputtering was carried out 
until no traces of oxygen or carbon was detected. The 
atomic diffusion was driven by varying the temperature of 
the samples. The bilayer sample and alloys were heated 
to 500 °C at rate 30 deg/min for different times and then 
cooled close to room temperature. 

The photoemission spectra were collected using both 
conventional x-ray photoelectron spectroscopy (XPS) 
(PHI ESCA 5400 Electron Spectrometer, Perkin Elmer) 
with non-monochromatic A1 Kq radiation at home lab¬ 
oratory and synchrotron radiation excited hard X-ray 
photoelectron spectroscopy (HAXPES) with high kinetic 
energy (HIKE) experimental station^^ at KMC-1 beam¬ 
line at Helmholtz-Zentrum Berlin (HZB), Bessy H. VG 
Scienta R 4000 electron analyzer, modified for electron 
kinetic energies up to 10 keV and highresolution double¬ 
crystal monochromator were used. To obtain photon en¬ 
ergies from 2300 to 7300 eV, Si(lll), Si(311) and Si(422) 
crystals of the monocromator were selected for the pre¬ 
sented measurements. The X-ray incidence angle was 
approximately 4° in every experiment and the photo¬ 
electrons were detected in normal emission. The energy 
scale was calibrated using Eermi level of the samples 
and Au 4f spectra of the calibration sample. The ad¬ 
justable photon energy range from about 2 keV to 10 
keV makes it possible to study photoelectrons with high 
kinetic energy which increases their inelastic mean free 
path (IMFP) making the HAXPES technique bulk sensi¬ 
tive as a comparison to surface sensitive laboratory XPS 
or soft x-ray range synchrotron radiation. Thus bulk sen¬ 
sitive investigations of atomic concentrations and chem¬ 
ical state of compound elements are possible without al¬ 
tering the original sample structure or chemistry by sput¬ 
tering. Also depth-profiling can be done by measuring 
core-levels with different binding energies and thus pho- 


4 


toelectrons with different IMFPs or by exploiting differ¬ 
ent sampling depths of specific core electrons by adjusting 
the photon energy of radiation. The high photon energy 
range makes it possible to measure photoelectrons with 
very high binding energy (low kinetic energy), for exam¬ 
ple Cr or Fe Is at 5990 and 7110 eV which enables also 
surface sensitive studies to be performed using HAXPES. 

In addition we have used conventional Auger elec¬ 
tron spectroscopy (AES) with Physical Electronics Model 
Cylindrical Mirror Analyzer to carry out depth profiling 
thorough the whole Fe/Cr bilayer before and after heat¬ 
ing until Si substrate is reached. The depth profiles of 
the samples were obtained using Ar+ sputtering (3000 V, 
4x4 mm^ area, about 8 mPa Ar+ pressure and 2 fiA 
sample current). After every 1 minute sputtering cycle 
O KLL, Cr LMM and Fe LMM spectra from 400 to 760 
eV were measured. The O and Cr spectra are overlap¬ 
ping and to get reliable quantitative information we used 
reference spectra of Fe, Cr and O (AlO) measured us¬ 
ing same equipment and parameters. The AES profiles 
were obtained from differentiated spectra of Fe/Cr bilay¬ 
ers and reference samples after background subtraction 
using factor analysis delivered by Principal Component 
Analysis (PCA) of Casa XPS 2.13.16 (Fig. 10). 

III. RESULTS AND DISCUSSION 

We begin the structural analysis of Fe-Cr alloys by us¬ 
ing the ab initio data to formulate some basic atomic 
models of the surface segregation and bulk precipitation 
of Cr. Next step is to perform MC simulations for larger 
systems to obtain results for more realistic cases. Finally, 
we use photoelectron spectroscopy to experimentally de¬ 
termine the segregation and precipitation in Fe/Cr dou¬ 
ble layer and Fe-Cr alloys as a function of annealing time 
and concentration. 


A. Initial relaxation in the surface region 

Besides on external conditions, the structural evolution 
of an initially homogeneous Fe-Cr alloy depends crucially 
on the concentration of the alloy and the relative mag¬ 
nitudes of the atomic diffusion rates at different spatial 
regions of the alloy. Using the ab initio data shown in 
Table I and in Figure 1 we can make predictions for the 
structural evolution of initially homogeneous Fei-xCi'a, 
alloys (same concentration throughout the whole bulk 
and surface regions). Since atomic diffusion rates are 
usually significantly higher in the near-surface regions 
than in the bulk^^ it is natural to split the consider¬ 
ation into initial surface relaxation and more retarded 
bulk relaxation. For non-homogeneous alloys, the struc¬ 
tural evolution depends also on the gradients of the Cr 
concentration. 

We begin our analysis by considering the relaxation 
of the near-surface region connected to the bulk reser- 
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TABLE I. Calculated (EMTO) bulk mixing enthalpy (Emix) 
and the bulk to surface (surface layer: si, second layer: S 2 ) 
segregation energy of Cr ) of homogeneous Fe-Cr 

as a function of atomic % of Cr in bulk (cb). Cb = 3 at.% cor¬ 
responds to the minimum of Emix. Negative (positive) mix¬ 
ing enthalpy means stable (metastable) bulk phase; negative 
(positive) E£g/’“^° means a driving force to a Cr containing 
(pure Fe) surface. The two upper rows are calculated using 
8 atomic layer slab in our previous work^'^ and the two lower 
rows are the present results of 12 atomic layer calculations. 



FIG. 1. (Color online) Calculated (EMTO) effective chem¬ 
ical potential A/r = /rpe — fJ-Cr of the surface atomic layer 
(red squares), second atomic layer (blue diamonds), and bulk 
(black circles) as a function of Cr content of the homogeneous 
alloy. The open symbols are calculated points, the curves are 
guides to the eye. In the present calculations we have used 
larger unit cells, 12 atomic (100) layers separated by vacuum 
of thickness equivalent to 6 atomic layers, compared to the 
previous work^”* where 8 atomic layers of metal slab was sep¬ 
arated by vacuum of thickness equivalent to 4 atomic layers. 


voir with fixed concentration. According to the ah initio 
segregation energy Si and 

S 2 refer to the surface and second atomic layer, respec¬ 
tively, Fig. 1), the initial driving force < 0) 

within the two surface atomic layers pushes Cr atoms 
from the second layer to the surface layer. This driving 
force is increased by a factor of three, from ^ 98 meV 
to ^ 295 meV, when the Cr concentration of the alloy 
increases from 0 at.% to 10 at.%. As Fig. 1 shows, the 
driving force on Cr atoms between the bulk and the sur¬ 
face atomic layer is for low Cr alloys from the surface to 
the bulk, but when the Cr concentration exceeds ^ 10 
at.% this driving force turns in the opposite direction 
pushing Cr atoms from the bulk to the surface. For all 
investigated alloy concentrations the second atomic layer 
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forms a diffusion barrier for a Cr atom to move from 
the bulk to the surface atomic layer. The barrier starts 
from '^310 meV at 0 at.% Cr and decreases to ~ 250 
meV at 10 at.% Cr. This can be compared with other 
ab initio calculations^^ using 36 atom supercell includ¬ 
ing one Cr impurity atom and giving 355 meV barrier 
for the second atomic layer. From Fig. 3 we see also 
that the model consisting of homogeneous bulk and two 
surface atomic layers describes the energetics of Fe-Cr 
(100) surfaces quite well because the segregation energy 
is practically converged to the bulk value at the third 
atomic layer. 

To sum up, our ab initio calculations predict a Cr de¬ 
pleted surface atomic layer for Fe-Cr alloys below 10 at.% 
Cr, a stable Cr enriched surface atomic layer within the 
bulk Cr content between 10 and 18 at.% (a typical sit¬ 
uation in many commercial steel grades), and a Cr con¬ 
taining stable surface beyond 18 at.% Cr. The second 
atomic layer is predicted to be depleted in Cr. Since 

increases with increasing Cr surface concen¬ 
tration, there exists a concentration dependent upper 
limit for the Cr content at the surface, as posed by the 
condition = 0.^^ 

B. Bulk-surface relaxation 

Considering the thermodynamic ground state of the 
whole bulk-surface Fe-Cr system, the bulk part should 
be relaxed too. The mixing enthalpy of Fei_a;Cra; (Ta¬ 
ble I) suggests that if the Cr concentration of the alloy 
exceeds ^10 at.%, the bulk part of the alloy has a ten¬ 
dency to transform to an a-a’ phase separated system: 
Cr-rich precipitates immersed in the Feo.gyCro.os alloy. 
The Cr-rich a’ precipitates are expected to avoid the con¬ 
tact with the surface because the surface energy of Cr is 
higher than that of Fe. Therefore, a low-Cr zone under 
the surface is expected to be formed driving the initially 
formed Cr-containing surface back to the pure Fe surface. 
However, this can happen only in vacuum. In ambient 
conditions the surface is expected to oxidize rapidly. Be¬ 
cause the Cr affinity to oxygen is much higher than that 
of Fe, CrgOs islands are expected to be formed on the 
surface. Due to the driving force of Cr to enrich to the 
metal/CrgOa interface^® these islands can grow until the 
uniform protective oxide layer is formed on the surface. 
Therefore, Cr at the surface is bound to an oxide form, 
and the surface of Fe-Cr is practically in an inert state 
during the retarded a-a’ phase separation in the bulk. 

C. Monte Carlo simulations 

The predictions of the basic properties of Fe-Cr sys¬ 
tems by the potential model are shown in Figs. 2, 3, and 
4 where the surface and segregation energies are shown. 

The energies of the (100) and (110) surfaces calculated 
by the 2BEAM are considerably lower when compared 


with the ab initio results. However, the differences be¬ 
tween the energies of iron and chromium are similar when 
comparing the 2BEAM and ab initio results. Further¬ 
more, all calculated energies differ from the experimen¬ 
tal values, which are also scattered and based partly on 
semiempirical estimates. 

Surface segregation energies for different alloy concen¬ 
trations were calculated as averages of 1000 random alloy 
samples. Segregation energy was defined as the energy 
difference between configurations where the Cr atom was 
in the center of the simulation box and when the atom 
was in one of the near surface atomic layers. One should 
note that the variation of the segregation energy in dif¬ 
ferent samples was large (in the range 0.1-0.2 eV) com¬ 
pared to the energy itself. The error bars (barely visible) 
in Figs. 3 and 4 are the errors of the mean. The figures 
show that there is a barrier for chromium atoms to seg¬ 
regate to alloy surface. However, for alloys of chromium 
concentrations in the range of 5-20 at.% the segregation 
energy to the second atomic layer is negative predict¬ 
ing Cr segregation to the second layer. Comparing the 
2BEAM results with the ab initio data of Levesque et 
al.^^ at the Cr impurity level the 2BEAM model predicts 
larger (smaller) segregation energy for the surface layer 
(second layer). On the other hand, the EMTO segrega¬ 
tion energy (Table I) is close to the 2BEAM value for 
the surface layer but close to the result of Levesque et 
al.^^ for the second layer. Fig. 4 shows essentially the 
same data as Fig. 3 but plotted along the Cr concentra¬ 
tion axis. Here the ab initio data is from Table 1. As 
Fig. 4 shows the EMTO results predict Cr segregation 
to the surface when Cr concentration exceeds 10 at.%, 
whereas the MCMD results predict that the segregation 
of Cr to surface is prevented by an energy barrier but 
Cr segregation is expected to the second atomic layer 
within the range of 5-20 at.% Cr concentration. One 
should also note that there is a discrepancy between the 
two ab initio results for surface layer of pure iron: 0.078 
eV (Ref. 21) vs. 0.216 eV (Table I). One notices from 
Figs. 3 and 4 that for the segregation energy Egegr the ah 
initio methods predict stronger oscillations as a function 
of layer position or the concentration than the MCMD 
method. This can be partly related to the difference of 
model systems used in the calculations. In the ab initio 
calculations there are only one type surroundings for a Cr 
atom in a specific layer whereas in MCMD calculations, 
in principle, every Cr atom has a different surrounding. 
Therefore the MCMD results are ’averaged’ which possi¬ 
bly leads to reduced oscillations in E^egr- 

The difference between the results obtained by various 
computational approaches is expected to be mainly due 
to the specific approximations and implementations of 
the computational methods. The EMTO results were ob¬ 
tained relaxing the volume of the alloy uniformly (with¬ 
out local atomic relaxations) and random occupation of 
atomic sites was simulated using the CPA, i.e. the model 
system in our case is a homogeneous alloy. Levesque 
et al.^^ used supercell approximation with fully relaxed 
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atomic coordinates to model an impurity Cr atom. Their 
model system consisted of one Cr atom in the unit cell. 
The average overall bulk concentration in their calcula¬ 
tions corresponds few at.% Cr and the planar concen¬ 
tration of Cr impurity is 25 at.% for (2a 2a 4a) slab (a 
being the lattice constant). For layer relaxations of pure 
Fe they reported —0.02 A for the surface layer and 0.04 
A for the second layer. A recent investigation based on 
the same ab initio method^^’^® found —0.002 A for the 
surface and 0.04 A for the second layer relaxations. This 
indicates that although there is a good agreement for the 
subsurface layer the top layer results are very different in 
these two pseudopotential ab initio investigations. Here 
we report similar behaviour for the segregation energies 
too. The differences between the results of Punkkinen 
et al.^"^’^® and Levesque et al.^^ could be partly due to 
the smaller slab and vacuum thicknesses and the using 
of constant-volume atomic relaxation in the latter inves¬ 
tigation. The semiempirical potentials in the MCMD 
method are optimized mainly for bulk properties, which 
somewhat limits its feasibility for surface studies. 

Since the accurate modeling of the surfaces and large 
bulk systems of Fe-Cr alloys at different temperatures 
is beyond the reach of any single computational tech¬ 
nique we are forced to use several theoretical approaches 
to get the comprehensive overall picture. The results of 
Levesque et al.^^ are considered to describe the impurity 
in bulk, but due to the relatively small cell size these data 
cannot account for the impurity effects on the surface. 
Considering surface regions at higher Cr concentrations 
the EMTO method is expected to perform better and 
for large bulk systems, as a function of temperature, the 
MCMD is the natural choice. 



Surface 

FIG. 2. (Color online) Surface energies of low index surfaces 
of bcc iron and chromium as predicted by the 2BEAM model. 
Results of ab initio calculations (filled squares) are from Ref. 
34 and 35. Experimental values (open triangles) are from Ref. 
36. 

Fig. 5 shows the near-surface concentrations from the 
surface layer (layer 1) up to the fifth atomic layer (layer 
5) as a function of bulk chromium concentration Cbuik- 



1 2 3 4 5 

Layer index 


FIG. 3. (Color online) Segregation energy of a chromium 
atom in iron chromium alloy in the near surface atomic lay¬ 
ers for different chromium concentrations as predicted by the 
2BEAM model (solid lines) and ab initio (dashed lines) cal¬ 
culations (Table I). Curves are shifted by the marked amount 
for better visibility. The ah initio result for the pure iron from 
Ref. 21 is plotted with dotted line and open symbols. 


At 300 K the surface layer is exclusively occupied by 
Fe atoms up to ^5 at.% Cr in bulk. At that point the 
Cr concentration at the surface jumps slightly upwards 
(to ^ 0.3 at.%). With increasing temperature the qual¬ 
itative shape of the concentration curve of layer 1 does 
not change appreciably, but the height of the jump gets 
larger, however, staying considerably lower than the av¬ 
erage Cr concentration in bulk (thin gray line in Fig. 5). 
The MCMD result for the bulk concentration threshold 
of the Cr containing surfaces (cbuik'^S-G at.%) compares 
reasonably well with the ab initio results (8-9 at.%).^^ 
However, the ab initio investigations predict the concen¬ 
tration of Cr at the surface to exceed the bulk value 
whereas in MCMD simulations the Cr concentration at 
the surface stays below the bulk value. This is related 
to the fact that the 2BEAM model predicts a strongly 
positive segregation energy for a Cr atom in the surface 
layer. 

The plateau beyond ^^10 at.% Cr in bulk seen in the 
concentration curves of the layer 1 suggests that the Cr 
content in nearby atomic layers have reached a certain 
saturation value. This is what actually happens as can 
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FIG. 4. (Color online) Segregation energy of the three near 
surface atomic layers as a function of bulk Cr concentration 
as predicted by the 2BEAM model. The ah initio data is the 
one presented in Table I. 

be seen in the panel of layer 2: the Cr content satu¬ 
rates to the ^^13 at.% value at the bulk concentration 
of Cbuik~8 at.%. With increasing depth from the sur¬ 
face the layer resolved concentration curves gradually ap¬ 
proach the average bulk concentration line, as expected. 
However, one should remember that the profiles shown in 
Fig. 5 are averaged concentrations parallel to the surface 
plane direction. The bulk part of Fe-Cr is expected to 
be a-a' phase separated at higher bulk concentrations, 
where the thermodynamically optimal state consists of 
Cr-rich precipitates (ccr = 80 — 90 at.%) in a homoge¬ 
neous Fe-Cr alloy containing few per cents of Cr. Since 
the surface energy of Fe-Cr is in this model minimized 
by an Fe surface, the probability of finding a Cr rich 
precipitation in the near-surface region is low. This is 
also clearly seen in Fig. 5. The Cr concentrations of the 
third and fourth layers level off to values below the av¬ 
erage bulk value when the Cr content in bulk is above 
the a-a’ phase separation threshold (average Cr content 
in bulk 6, 10, and 15 at.%, for the 300, 500, and 700 K 
simulations, respectively). In the case of layer 2 the in¬ 
crease of the concentration above the average bulk value 
at Cbuik = 3 — 13% can also be attributed to the negative 
segregation energy of a Cr atom in the second layer at 
these bulk concentrations. 

In Fig. 6 the cross sections of the MCMD simulation 
box are shown for concentrations Cb = 5, 10, 15, and 
20 at.% at simulation temperatures 300 K and 700 K. 



0 2 4 6 8 10 12 14 16 18 20 

^bulk 

FIG. 5. (Color online) Chromium concentration of near¬ 
surface atomic layers as a function of the average bulk con¬ 
centration. Gray lines show the one-to-one relation between 
the bulk and layer concentrations. Note the different y axis 
scale in the case of layer 1. 


The low probability of Cr rich a’ precipitates touch¬ 
ing the surface is clearly demonstrated. Adding Cr in¬ 
creases the number and size of Cr-rich precipitates leav¬ 
ing the rest of the alloy in a homogeneous low-Cr state 
in agreement with recent theoretical and experimental 
investigations.Fig. 7 shows the concentration pro¬ 
files corresponding to the MCMD simulations shown in 
Fig. 6. One can see that when the bulk concentration 
Cbuik is below ^10 at.% the surface depth profile of the 
Cr concentration shows only the short-period screening 
effect due to the surface perturbation leveling off to the 
bulk concentration in deeper layers. For concentrations 
Cbuik ^ 10 at.% long-period oscillations due to the Cr- 
rich precipitates show up deeper in the sample. An al¬ 
most fixed ^13 at.% Cr value in the second atomic layer 
is observed due to the interactions of Cr-rich precipitates 
and the surface in line with the ab initio predictions. 
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FIG. 6. (Color online) Snaphots of the simulated 
Fei_a;Cra; systems at different temperatures and with differ¬ 
ent chromium concentrations x. Figure shows slices of thick¬ 
ness 15 A in the (110) direction. Only chromium atoms are 
shown with color-coding according to potential energy (blue 
colour representing low potential energy). The top and bot¬ 
tom facets of the simulation box are open surfaces. 


Fig. 8 shows similar depth profiles as Fig. 7 but with 
the Cr precipitates removed. This was done by discard¬ 
ing all the atoms that had the average atomic type index 
smaller than 1.5 where 1.0 corresponds to pure Fe and 
2.0 to pure Cr. We see that in the bulk part the con¬ 
centration reaches the solubility limit of the Cr in Fe (6, 
10, and 15 at.%, for the 300, 500, and 700 K simulations, 
respectively). Concentrations of the surface layers are 
essentially the same as in Fig. 7 because the chromium 
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FIG. 7. (Golor online) Chromium concentration depth profiles 
for different average chromium concentrations and simulation 
temperatures. Curves for 0.1 % are multiplied by 10 for better 
visibility. 


precipitates do not extend to the surface layers. With 
the effect of the precipitates removed. Fig. 8 also shows 
clearly the range and the decay rate of the concentration 
fluctuations due to the screening of the surface perturba¬ 
tion. The low Cr content in the first layer is overcompen¬ 
sated by the layer 2, which is back compensated by the 
layers 3 and 4 etc. As expected, the compensation ampli¬ 
tude decreases with increasing distance from the surface 
and the concentration fluctuations are less pronounced 
at higher temperatures. 

Considering the experimental investigations of segrega¬ 
tion and precipitation in Fe-Cr it would be beneficial to 
study, in addition to alloys, also Fe/Cr layer structure. In 
the layer structure one can study the evolution of Fe-Cr 
as a function of Cr concentration in a more transparent 
way because the density gradient is a one-dimensional 
function. Similarly in this case the segregation of Cr in 
the surface of the Fe layer can also be investigated in 
a controlled way. These benefits in mind, we finish our 
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FIG. 8. (Color online) As in Fig. 7 bnt with chrominm pre¬ 
cipitates removed. Removing was done by only taking into 
account atoms whose average type during the last 60000 sim¬ 
ulation steps was smaller than 1.5 where value 1.0 corresponds 
to pure Fe and 2.0 to pure Cr. 


simulations by considering Fe/Cr layer structure. Fig. 9 
shows the results of the simulations of the Fe/Cr inter¬ 
face. The dashed curves are the initial states in which 
there exist pure Cr and Fe sections in the simulation 
cell, whereas the solid lines show the situation after the 
simulation. The results show that the thermodynamical 
ground state of the investigated system consists of two 
different parts: low Cr section where the average Cr con¬ 
centration is ~ 6 at.% at 300 K 12 at.% at 700 K, 
which compares well with experimental results of about 
14 at.% at 773 K®) and the excess of the Cr remains in the 
original Cr part. In the case of temperature of 700 K and 
original concentration of 10 at.% the whole Cr layer is dis¬ 
solved into the Fe. This is expected as the solubility of 
Cr in Fe at this temperature is over 10 at.%. At the lower 
temperature of 300 K the original layer is partly dissolved 
leaving precipitates with mainly (110) facets. Also, for 
the 700 K case with 20 at.% chromium the layer structure 


changes into a single precipitate again with (110) facets. 
The appearance of (110) facets is in line with ab initio 
interface calculations®’’ which show that (110) interface 
has the lowest energy which compares well with the sur¬ 
face energies in Fig. 2. The 300 K, 50 at.% case shows 
oscillations of the concentration profile which are caused 
by the similar compensation effect as on the surface and 
partly by the faceting of the interface. These oscillations 
are washed out when temperature is raised to 700 K and 
the positions of the atoms become more random. 
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FIG. 9. (Color online) Gr concentration profiles for the layer 
widths corresponding different total Cr concentrations simu¬ 
lated in different temperatures. Thin dashed lines show the 
initial Cr profile. 


D. Segregation, precipitation and oxidation in 
Fe-Cr systems 

1. Cr segregation in Fe/Cr bilayer and Fe-Cr alloys 

Experimentally the segregation of Cr can be demon¬ 
strated by comparing the concentrations of Cr and Fe 
atoms as a function of the probing depth. We have used 
two different methods to carry out depth profiling. In 
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AES the topmost atoms are removed by ion gun layer 
by layer and after every sputtering cycle the Cr and Fe 
Auger electron spectra (Fig. 10) are collected and the in¬ 
tensities of the compound elements are compared. How¬ 
ever, due to ion bombardment sputtering is considered as 
a destructive method and it may have unwanted effects 
on some concentrations values. For example, preferen¬ 
tial sputtering or reduction of oxides can be a drawback 
in some studies. Preferential sputtering is not a prob¬ 
lem in case of Fe and Cr but the analysis of oxide layer 
can be. This technique is not ideal for chemical analysis 
which is why HAXPES spectra measured using different 
photon energies (different IMFP of photoelectrons) was 
important to collect as well. Due to the limited time 
for synchrotron radiation experiments it is impossible 
to measure HAXPES spectra layer by layer but use of 
three selected photon energies and analysis of different 
core-level spectra bring us information on the concentra¬ 
tions of Fe and Cr at different depths from the surface. 
The investigated Fe/Cr double layer and Fe-Cr alloys 
are known to have polycrystalline structure containing 
grain boundaries and compound intermixing and diffu¬ 
sion around Fe/Cr interface and oxygen diffusion into the 
layers^® which can complicate the investigation of Cr seg¬ 
regation and interpretation of the results. AES together 
with XPS was used to carry out preliminary diffusion and 
depth profile experiments. AES depth profiles that were 



FIG. 10. (Color online) Example of an AES spectrum of 
Fe/Cr interface (o curve) of Fe/Cr/Si sample. The sum spec¬ 
trum (dashed curve) was obtained by Principal Component 
Analysis (PCA) of Casa XPS 2.13.16.®® 

measured from untreated and heated Fe/Cr/Si layer sam¬ 
ples are presented in Fig. 11. Annealing was performed 
in XPS located in the same room than AES but missing 
the possibility to transfer the sample in UHV between the 
spectrometers. XPS spectra were measured before and 
after every heating period (Fig. 11) which were carried 




Depth from the surface [nm] 



FIG. 11. (Color online) Left hand side presents the Fe 2p, Cr 
2p and O Is XPS spectra of a Fe/Cr bilayer sample before 
(topmost spectrum) and after annealing. These spectra were 
measured using A1 Ka twin anode. Photoemission spectra 
clearly demonstrate the increase of Cr 2p intensity as a func¬ 
tion of annealing time confirming the Cr diffusion into the Fe 
layer. The upper AES profile on the right is taken over the 
Fe/Cr interface before any heating was carried out starting 
sputtering from the surface represented by upper most XPS 
spectrum (blue) on the left. The lower AES depth profile on 
the right was taken after all the heating treatments starting 
from the surface concentration presented in the lowest (red) 
XPS spectrum on the left. 

out at 500 °C for 5, 15 and 40 minutes giving total heat¬ 
ing time of 60 minutes. The pressure of the preparation 
chamber where the heating was performed was 3 x 10“® 
Torr at maximum. Sample was cooled close to room tem¬ 
perature before measuring the spectra. Annealed sample 
was transferred to AES through air and two minutes of 
sputtering was performed before starting the depth pro¬ 
filing to get rid of the carbon contamination. The upper 
depth profile in Fig. 11 nicely presents the well separated 
Fe and Cr layers. It must be noted that the signal in the 
topmost XPS spectra at the same binding energy as Cr 2s 
would appear (around 696 eV) is entirely due to the non- 
monochromatic radiation used in laboratory XPS. About 
10 % oxygen concentration in the Fe layer is comparable 
to the topmost XPS spectrum in Fig. 11. The content 
of oxygen through the whole sample is quite high but 
can be explained by the fact that polycrystalline samples 
contain number of grain boundaries or dislocations which 
can enhance the oxygen diffusion or diffusional transport 
of Cr towards the surface.Also the pressure during the 
sample growth and sputtering with Ar’*' ions can have ef¬ 
fect on the oxygen concentration within the layers. The 
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oxygen content increases in phase with Cr content in both 
profiles which in addition to the above mentioned reasons 
is related to higher affinity of Cr for oxygen and Cr react¬ 
ing with inward diffused oxygen.After heating the 
Fe/Cr bilayer for 60 min at 500 °C the layered structure 
has been destroyed and the Cr concentration within the 
first 40 nm is about 30 % increasing after that up to 
60 % before any signal from Si substrate was observed. 
The oxygen concentration for the annealed sample is less 
than 5 % until fast decrease in Fe and slow increase of Cr 
concentration takes place at approximately same depth 
from the surface where the original Fe/Cr interface was 
detected before heating. Oxygen concentration increases 
about the same percentage value as Cr concentration but 
more rapidly. This is expected since the diffusivity of O 
in bcc Fe^^ is much higher than diffusivity of Cr in Fe.^® 
Most probably the reason for low oxygen concentration 
within the first 40 nm in annealed double layer sample 
is due to the formation of protective passive layer on 
the surface because of increased Cr concentration. The 
concentration at the uppermost atom layers of annealed 
sample is not visible in Fig. 11 since some sputtering was 
done prior to the presented data to remove the contami¬ 
nation (carbon) that was absorbed on the surface during 
the sample transfer from XPS to AES. 

To investigate the Cr segregation and also precipitation 
in Fe-Cr systems and to test our theoretical predictions 
HAXPES technique was used to monitor the concentra¬ 
tion profiles of the Fe/Cr double layer and chemical state 
of Fe and Cr atoms in two Fe-Cr alloys as a function 
of photon energy and heating time. Investigations were 
performed with photon energies from 2300 to 7300 eV in 
order to study the concentration profiles of the chemical 
components as a function of the probing depth. Prior to 
HAXPES experiments several annealing tests were per¬ 
formed with laboratory XPS to find the right annealing 
temperatures and time. The ideal temperature for Cr 
segregation in Fe/Cr/Si turned out to be 500 °C. At this 
temperature Cr diffused towards the surface in a few min¬ 
utes time so that we were able to detect it with XPS. In 
the proper HAXPES segregation experiments the bilayer 
sample was heated in several steps as described above to 
gradually follow the Cr diffusion into the Fe layer (Fig. 
12). Before starting the HAXPES segregation study the 
Fe layers of Fe/Cr/Si were cleaned in the UHV of an¬ 
alyzer chamber of HIKE spectrometer by annealing at 
150 °C followed by 20 min Ar+ sputtering. To start the 
Cr segregation investigations the quality of the Fe/Cr 
double layer structure, i.e. the absence of Cr in the Fe 
layer, was checked by scanning the binding energy range 
of Cr 2p and Cr Is core levels before any heating was per¬ 
formed. Initially no Cr was found using 2300 eV, 4000 
eV and 7300 eV photon energies i.e. experimental condi¬ 
tions that correspond to IMFP (A) values from 2 to 7 nm 
and thus sampling depths (3A) of approximately 9-20 nm 
depending on the oxidation and core-level studied.It 
is worth mentioning that no Si was found in the spectra, 
either. This indicates the proper double layer structure 



FIG. 12. (Color online) Overview HAXPES spectra of the 
Fe/Cr bilayer sample measured with 2300 eV photon energy 
before and after two heating treatments. In addition to the Cr 
2s and Cr 2p spectra appearing after annealing, the increased 
intensity of Cr 2p as a comparison to Fe 2p after different 
annealing times demonstrates the Cr diffusion towards the 
surface. The intensity of Fe 2p3/2 of the presented spectra is 
normalized to one after setting background to zero. 


of the sample, complete Fe layer on top of a Cr layer. 
Due to sputtering the thickness of Fe layer before any 
heating treatments was less than in the AES studies be¬ 
ing still more than 20 nm. After the cleaning procedure 
sample was heated to 500 °C sequentially for 5, 15 and 40 
minutes, total heating time being 5, 20 and 60 minutes, 
respectively. The spectra were measured after each heat¬ 
ing period when the sample was cooled close to room 
temperature. The gradual changes of the atomic con¬ 
centration profiles were observed with heating time (Fig. 
12). These spectra nicely demonstrate the appearance 
of Cr spectra for heated samples. Before heating no sig¬ 
nal at Cr 2p binding energy was observed but already 
after 5 min heating Cr 2s and Cr 2p signals appeared 
as shown in Fig. 12. The cross sections or spectrome¬ 
ter transmission are not taken into account in Fig. 12. 
which slightly underestimates the real intensity of Cr 2p 
and 2s photoemission signals relative to that of Fe 2p. 
The bulk concentrations of O, Cr and Fe after short and 
longer heat treatments are shown in Fig. 13. Intensities 
presented in Fig. 13 are calculated by comparing the area 
of O Is, Fe 2p and Cr 2p spectra where Scofield cross- 
sections,^® spectrometer transmission and other experi¬ 
mental parameters have been taken into account. The Cr 
concentrations presented in Fig. 13 are in line with the 
AES depth profile in Fig. 11. The oxygen concentration 
for heated sample derived from the HAXPES spectra is 
higher than the AES profiles give but this is most prob¬ 
ably due to the pressure difference during the heating 
and measurements. Measurement with 2300 eV photon 
energy gives the average Cr concentration within 5 or 6 
nm thick layer to be close to 15 at. % (Fig. 13). The 
concentrations in the 4000 eV photon energy measure¬ 
ments are similar but when compared to the more bulk 
sensitive, 7300 eV photon energy case, the Cr concen¬ 
tration is much lower than the Fe concentration. The 
increasing Cr concentration in the direction from bulk to 







12 


the surface within the Fe layer of the heated sample is 
in line with our earlier calculations.In Fe0.g5Cr0.05 
and Feo.85Cro.15 alloys the Cr segregation was not so ev¬ 
ident. In bulk (comparison of Fe 2p and Cr 2p spectra) 
almost no change in Fe/Cr intensity ratio was detected 
but more surface sensitive measurements of Cr and Fe Is 
photoelectrons show some increase in Cr intensity as a 
result of heating, especially in the oxidized Cr. 


Fe/Cr ratio in Fe/Cr/Si 

■ Cr[%] BFera aora 
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FIG. 13. (Color online) Concentrations of Cr, Fe and O in 
Fe/Cr/Si. Relative intensities are calcnlated by analyzing the 
area of Cr 2p, Fe 2p and O Is spectra. 


2. Oxidation 

Even after 5 min annealing at 500 °C at pressure of 
about 5 X 10“® Torr the Fe/Cr bilayer sample and Fe-Cr 
alloys were evidently oxidized enabling the investigation 
of the ratio of alloy Cr and oxidized Cr. Oxidized Cr is 
well resolved on the high binding energy side of the alloy 
Cr (Fig. 14). A higher amount of oxidized Cr in Fe/Cr bi¬ 
layer was observed in the more surface sensitive (photon 
energy=2300 eV) Cr 2p HAXPES experiment (Fig. 14a) 
compared to the 4000 eV or 7300 eV measurements^® 
indicating increasing Cr concentration towards the open 
surface also in the initial oxidation experiments. The 
spectra show also that, especially at room temperature 
range, oxygen penetration deeper into the metal is slow. 

Fig. 14 presents Cr 2p photoemission spectra of bi¬ 
layer sample after 5 and 60 minutes heating at 500 °C 
measured with 2300 eV photon energy (a) together with 
Cr 2p spectra of Fe0.95Cr0.05 and Feo.85Cro.15 alloys be¬ 
fore and after heating 10 minutes at 500 °C (b and c, 
respectively). The spectra that were measured after 5 
min annealing of Fe/Cr/Si (a) and before annealing of 
alloys (b, c) present a fit with three components. The 
two components, Cr and Cr alloy can be resolved due to 
the chemical shift between Cr atom surrounded mostly 
by other Cr atoms or Fe atoms. Thus segregation of Cr 





FIG. 14. (Color online) Cr 2p spectra of Fe/Cr bilayer sample 
after 5 (open circles) and 60 min (solid line) heating (a) and 
Cr 2p spectra of Fe0.95Cr0.05 (b) and Feo. 85 Cro .15 (c) alloys 
before and after heating 10 minutes at 500 °C. The lower 
spectra in a, b and c are fitted using three components to 
represent signals from bulk like Cr, Cr from Fe-Cr alloy and 
oxidized Cr. 


and formation of Cr dominating areas can be estimated 
by analyzing the fitted spectra. The broad feature from 
approximately 575 to 580 eV in all the spectra is the 
spectrum of measured Cr203 reference sample detected 
with same parameters than the actual samples. The ap¬ 
pearance of this pronounced feature on the high binding 
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energy side of the bulk and alloy Cr makes it possible to 
distinguish the oxidized Cr froirr the unoxidized one. 

In Cr oxides the interaction between valence band 3d 
electrons and 2p vacancy core levels can create a num¬ 
ber of hnal states which is called multiplet splitting®°’®^. 
For example, Biesinger et al.®° have presented a fit for 
asymmetric Cr 2p photoemission spectra of Cr(III) ox¬ 
ide that consists of five components. In this study we 
have used Cr 203 reference spectrum to replace the mul¬ 
tiple feature fit that otherwise would have been needed to 
properly describe the complicated multiplet splitting of 
Cr(III) oxide®°. However, the Cr 2p photoemission spec¬ 
tra can’t be fitted satisfactorily by using only Cr 203 ref¬ 
erence spectrum and one asymmetric Voigt-line shaped 
feature to describe the unoxidized Fe-Cr alloy. An ad¬ 
ditional component is needed between Cr 203 and Cr al¬ 
loy features. This structure at approximately 574.4 eV 
binding energy originates from bulk like chromium (Cr- 
Cr) (labelled Cr in Fig. 14). The shift between the Cr 
alloy and Cr components in the fitted spectra is approxi¬ 
mately 0.9 — 1 eV, which is close to the value we observed 
in our HAXPES measurements of pure, unoxidized Cr 
and Cr 203 reference samples. The asymmetry parame¬ 
ter and full width at half maximum (FWHM) value used 
for Cr and Cr alloy line shape in the fitting procedure 
were derived from the HAXPES spectra of Cr bulk ref¬ 
erence. In addition to the experimental arguments first 
principles calculations within density functional theory 
using the complete screening picture®^ were performed to 
estimate the binding energy shift between bulk Cr and 
Fei-xCra; alloys of different concentrations. These cal¬ 
culations gave a negative binding energy shift for Cr 2p 
core level spectra which was between —0.1 and —0.5 eV 
for Cr concentrations x = 0.1 — 0.9 being largest in case 
of Feo, 8 Cro, 2 . 

Even though the samples were measured in ultrahigh 
vacuum the Cr oxidation was evident. However, almost 
no sign of Fe oxide was seen in the Fe photoemission 
spectra. As presented already in Fig. 14 intensity of Cr 
in Fe/Cr/Si further increased between the 5 and 60 min¬ 
utes heating and in the latter spectrum larger number of 
Cr atoms is detected as an oxide than in metallic form 
(Fig. 14a). During the 60 min annealing the amount of 
Fe decreases (Fig. 13) to lower level due to Cr segregation 
(Fig. 14a) and the oxidized form of Fe is not resolved even 
on the surface. As a comparison to unoxidized (sputtered 
until no traces of O or C were detected) bulk Fe refer¬ 
ence sample only the surface sensitive Fe Is spectrum 
of Fe/Cr/Si (not shown here) that was measured before 
any heating shows slightly increased intensity on the high 
binding energy side where oxidized Fe would show. 

The Cr 2p spectra of Fe 0 . 95 Cr 0.05 and Feo. 85 Cro .15 al¬ 
loys in Figs. 14b and 14c show how the concentration 
can affect the formation of oxide layer already at very 
low pressure (10“® Torr). These spectra were measured 
using 4 keV photon energy in which case the share of sur¬ 
face layer (topmost 5 nm) has here a smaller role than in 
the Fe/Cr/Si spectra of Fig. 14a. In the case of higher 


Cr bulk concentration alloy (Feo. 85 Cro. 15 ) most of the Cr 
atoms are oxidized even before any heating was carried 
out. The fraction of Cr 203 component of the total Cr 
2 p 3/2 signal area is 55 % for Feo. 85 Cro .15 and 23 % for 
Fe0.95Cr0.05 alloy. In Fe0.95Cr0.05 the heating seems to re¬ 
move part of the Cr 203 relative to unoxidized Cr making 
the Cr/Cr oxide signal ratio slightly increase. Still the 
amount of unoxidized Cr stays approximately the same 
despite the heating. Also here only the Cr atoms were 
oxidized. 


3. Cr solubility and precipitation 

While the total amount of Cr increases considerably in 
Fe/Cr/Si (Fig. 13) due to additional heating the share 
of oxidized Cr has the highest proportion. On the other 
hand, the Fe 2p spectra measured with photon energy of 
7300 eV reveal that with heating the Fe concentration 
in the bulk part of the Fe-layer changes only moderately 
compared to the surface region (Fig. 13), implying sim¬ 
ilarly moderate changes for the Cr content in the bulk 
part of the Fe-layer. As the MCMD simulations show, in 
the Fe/Cr double layer system chromium is dissolved in 
the iron layer up to the Cr solubility limit: 6 at. % at 300 
K and 12 at. % at 700 K (Fig. 8). Thus our MCMD result 
for the Cr solubility limit is in line with the present HAX¬ 
PES measurements of 5 and 60 minutes annealed Fe/Cr 
double layer showing that the Cr concentration is lowest 
in the Fe him between the underlying Cr layer and the 
surface region. This demonstrates the existence of the 
upper limit for the Cr solubility in Fe and therefore, in 
the present case, the Fe layer acts as a retarder for the 
Cr diffusion from the deeper lying Cr reservoir towards 
the surface and the Cr enrichment to the surface. 

The MCMD simulations of formation of Cr rich a' pre¬ 
cipitates at different temperatures with different Cr con¬ 
centration in Fei_a;Cra; systems are presented in Fig. 6. 
The experimental techniques used in this study are not 
suitable to study precipitation in such a resolution that 
the measurements could be carried out only for the areas 
where precipitates occur as a comparison to the alloy ar¬ 
eas. However, the binding energy of Cr atoms is very sen¬ 
sitive to their chemical state which is why combination of 
Cr bulk reference sample spectra and information based 
on hrst principles calculations using the complete screen¬ 
ing picture were used to fit the Cr 2p spectra (Fig. 14). 
The feature labeled Cr has higher binding energy than 
the component that describes the fraction of Cr atoms in 
Fe-Cr alloy. This difference depends on the Cr concen¬ 
tration in alloy being here approximately 0.9 eV referring 
to Cr concentration around 15 — 30 at.% which cause the 
largest negative binding energy shift as a comparison to 
100 % Cr. The feature labelled Cr originates from atoms 
that have more other Cr atoms than Fe (or O) atoms as 
nearest neighbors. In alloys where Cr concentration is 
as low as 5 or 15 at. % the origin of this Cr bulk like 
structure can be either formation of Cr rich layer under 
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the topmost surface atom layers due to segregation or 
precipitation of Cr. Cr can also have high concentration 
on the grain boundaries but then it would most likely be 
oxidized. According to the simulations in alloys with Cr 
bulk concentration exceeding 10 at. % (Fig. 6 ) Cr rich a’ 
precipitates can be found already at room temperature. 
In the experiments the relative intensity (area of the Cr 
component) increases considerably from Fe0.95Cr0.05 to 
Feo. 85 Cro .15 alloy being approximately 7 and 12 %, re¬ 
spectively. This trend suggests that at least part of the 
Cr signal can be caused by Cr segregated a’ phase, i.e. 
Cr precipitates in bulk, being in line with the theory. 

However, more detailed experiments and discussion of 
the bilayer sample and Fe-Cr alloys are still needed be¬ 
fore more thorough conclusions about precipitation and 
oxidation states of Cr after every heating step can be 
drawn. Here the fit is used to give a better picture how 
the method can be used to follow the rate of oxidation 
and progression of Cr segregation and to estimate the 
Cr concentration over the analyzed sample layer. When 
comparing the AES, XPS and HAXPES results besides 
resolution and probing depth of the techniques, it is es¬ 
sential to keep in mind the possible effect of sputtering 
on the chemical composition of the bilayers, the different 
rate of oxidation caused by differences in analyzer cham¬ 
ber pressures and possible treatment and transfer of the 
samples prior to the measurements. All these mentioned 
parameters were set to meet each other as well as pos¬ 
sible and the differences have been taken into account 
when analyzing the results. 


IV. CONCLUSIONS 

Due to the major challenges related to the investiga¬ 
tion of long-range properties of materials as a function of 
different internal and external parameters we have used 
in the present investigation a multi-scale and interdisci¬ 
plinary approach. EMTO method has been used to in¬ 
vestigate atomic-scale properties, MCMD method is ap¬ 
plied for exploring the large-scale bulk phenomena, and 
several spectroscopic techniques have been used to study 
properties related to kinetics and oxidation. 

The ab initio EMTO calculations predict that in the 
initially homogeneous Fei-^Cr^; the net driving force of 
the migration of Cr atoms is from the second atomic layer 
to the bulk for low Cr alloys (ccr ^ 10 at.%) and to the 
surface for moderate Cr alloys (ccr ^ 10 at.%). Com¬ 
paring the bulk and surface the calculations predict the 
driving force to be from the surface to the bulk (ccr ^ 10 
at.%) and from the bulk to the surface (ccr ^ 10 at.%), 
i.e. the ab initio simulations predict Cr containing sur¬ 
faces when Cr concentration exceeds ~10 at%. 

Monte Carlo molecular dynamics simulations predict 
that the formation of Cr-rich precipitates affect the ther¬ 
modynamics of the surface layers of Ee-Cr alloy. Increas¬ 
ing the Cr concentration the excess chromium goes to the 
precipitates. Due to the Cr surface and segregation ener¬ 


gies the precipitates do not reach the surface layer and, 
consequently, there is no large thermodynamic driving 
force to push Cr atoms to the surface. In ambient condi¬ 
tions oxidation of Cr changes the picture. Surface Cr is 
oxidised and the oxide-alloy interface forms a sink to the 
Cr of the alloy phase. One should, however, note that 
these simulations only give the thermodynamic equilib¬ 
rium states. Kinetics of the alloy microstructure and 
surface oxidation can be expected to have effect on the 
surface structure. Simulations of Fe/Cr layer structure 
show the temperature dependent upper limit of the sol¬ 
ubility of Cr in Fe ~ 6 at.% at 300 K and 12 at.% at 
700 K. 

Experimental investigations concentrated on Cr segre¬ 
gation and precipitation together with initial oxidation of 
Fe-Cr systems. The segregation of Cr was nicely demon¬ 
strated in the spectra by annealing Fe/Cr double layer 
and following the diffusion of Cr towards the surface. 
This was done using both AES and HAXPES measure¬ 
ments. Besides comparison of the intensity of Cr and Ee 
core level spectra in Ee/Cr double layer and Fe-Cr alloys 
more detailed information of chemical state of Cr atoms 
was derived from the deconvoluted Cr 2p spectra. The 
fitting procedure suggested that part of the Cr atoms 
have bulk Cr like structure which can be connected to 
the formation of Cr rich precipitates where the nearest 
neighbours of Cr atoms are other Cr atoms instead of 
Fe (or O) atoms. The share of this Cr signal was much 
higher for Feo. 85 Cro .15 than for Fe0.95Cr0.05 alloy which 
is in line with the MCMD simulations. Initial oxidation 
of Fe-Cr systems is provably very demanding to investi¬ 
gate because of very fast Cr oxide formation in freshly 
cleaned samples even at very low pressure, in the UHV 
of analyzer chamber of spectrometers. Cr was the only 
compound that was oxidized during and after annealing 
and the formation of Cr 203 was verified by using the 
reference spectra and fitting procedure. 

In order to further improve the understanding of the 
properties and phenomena related to Fe-Cr in different 
ambient conditions the development trends such as MC 
simulations using potentials which properly take account 
for the surface effects (e.g. enhanced magnetic moments), 
as well as more extensive DFT simulations using effec¬ 
tive chemical potentials in kinetics would be advanta¬ 
geous. The developing HAXPES technique will also pro¬ 
vide more selective research methods in future. 
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